Gutzwiller density functional theory for correlated electron systems 
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We develop a new density functional theory (DFT) and formalism for correlated electron systems by tak- 
ing as reference an interacting electron system that has a ground state wavefunction which obeys exactly the 
Gutzwiller approximation for all one particle operators. The solution of the many electron problem is mapped 
onto the self-consistent solution of a set of single particle Schrodinger equations analogous to standard DFT- 
LDA calculations. 

PACS numbers: 



Over the last several decades, first principles total energy 
calculations using density functional theory based on the lo- 
cal density approximation (LDA) or generalized gradient ap- 
proximation (GGA) have been well developed into a theoret- 
ical tool with strongpredictive capability for a large number 
of materialsJIlS[3n3l. However, there are important classes 
of materials involving strongly correlated electrons ranging 
from high Tc superconducting compounds and various other 
transition metal oxide materials to f-electron elements bear- 
ing materials where the current LDA/GGA approaches fail 
in fundamental ways. There have been intensive studies on 
new approaches to remedy the situation such as LDA+uJl, 
LDA-dynamical mean field theory (DMF^JHHEIl, self- 
interaction correction-local spin density (SIC-LSD) llOll . and 



hybrid f unctionals II 1 lL 1 1 211 . Although these approaches par- 



tially address the issues related to the strongly correlated elec- 
tron systems, a comprehensive and generally accepted pre- 
dictive theory with the quality of LDA for normal metal, al- 
loys, and compounds is still lacking for materials containing 
strongly correlated electrons. To address this problem, we 
propose in this paper a new density functional theory that goes 
beyond LDA through a self-consistent solution of the many- 
body ground state using the Gutzwiller approximation lfl3l for 
interacting electron systems within a first-principles frame- 
work. 

The Hohenberg-Kohn theorem states that the ground state 
energy of an electron system is a functional of the electron 
densityl 1]. Kohn and Shamj2l took this a step further by ex- 
pressing the energy of a real system in terms of the energy 
of a fictitious system of non-interacting electrons that has the 
same density and same kinetic energy as the real system. This 
led to a system of non-interacting electrons moving in an ef- 
fective potential that can be solved through the iterative solu- 
tion of a set of one-electron Schrodinger equations within the 
LDA for the exchange-correlation energy. The current first- 
principles density functional calculations are based on this set 
of effective one-electron Kohn-Sham equations. It should be 
noted that the Hohenberg-Kohn density functional theorem 
is true for any electron system including strongly correlated 
electron systems. The failure of the Kohn-Sham approach for 
the strongly correlated electron systems, in our view, is largely 
due to the choice of non-interacting electron as the reference 
system. By carefully choosing a reference system which in- 



cludes the most essential strong correlations, yet still can be 
cast into a set of one-electron Schrodinger equations through 
variational principles, a new density-functional formalism for 
treating strongly correlated electron systems can be derived 
and implemented following the spirit of Kohn and Sham[2]. 

In the new approach we propose here, instead of defining 
the kinetic energy functional to be the kinetic energy of a 
non-interacting electron gas with the same density, we will 
define the kinetic energy functional to have a simple ana- 
lytical form corresponding to the frequently used Gutzwiller 
approximation for a sy stem of electrons with on-site-only 



correlations! 13, 14, 15, 3, 17, 18]. The expectation of any 



one particle operator (e.g. the electron kinetic energy and 
the electron density) under the Gutzwiller approximation can 
be expressed in terms of a non-interacting one-particle den- 
sity matrix with renormalized wei ght d ue to adjustments from 
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strong correlation effects [13, 
show that, within the Gutzwiller approximation, the inter- 
acting many-electron problem can be mapped onto a non- 
interacting system with an effective potential. The exact 
Coulomb interactions can be included for a pre-determined set 
of localized configurations while the local density approxima- 
tion is used for all the remaining exchange-correlation contri- 
butions. 

The choice of including electron correlations using the 
Gutzwiller approximation is motivated by previous work 
where it has been shown to interpolate well between the two 
regimes of strong electron correlation (large U limit) and 
small electron correlation (small U limit). Another impor- 
tance of the Gutzwiller approach is the correct description 
of highly correlated states near the Fermi level. LDA can 
be view as an extension of Hartree-Fock theory into density 
functional theory, we view our present scheme as an extension 
of quantum chemical couple-cluster calculations into density 
functional theory. In our scheme, the variational parameters 
are the single particle electron wavefunctions for the localized 
and delocalized electrons and the occupancy of the various lo- 
calized configurations at each atom in the unit cell. Like LDA, 
the formulation is from first principles with all Coulomb and 
exchange interactions determined self-consistently. There are 
no adjustable parameters. 

According to the Hohenberg-Kohn density functional the- 
orem, the ground state energy of a multi-electron system is a 
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functional of the electron density p 
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Instead of taking as reference a non-interacting electron gas 
with the same density as the exact many-electron system, the 
density functional in our present theory is determined by tak- 
ing as reference an interacting electron system: 
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Exc [p] 



p(r)v{r,r')p(r')d 3 rd 3 r' (2) 



The reference system is chosen to have the same electron den- 
sity p (r) as the ground state of the exact multi-electron system 
and to have a ground state wavefunction I^g) which obeys 
exactly the Gutzwiller approximation for all one particle op- 
erators 



N 



(3) 



m— 1 



In the Gutzwiller approximation, for each one-particle oper- 
ator acting on |^g)> we can define a corresponding renor- 
malized operator Oq, acting on the underlying Hartree-like 
wavefunction |\I/o)used m generating \^g) m the Gutzwiller 
approach such that: 
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Therefore we have 



O 



X - - (*o 



#0 



*0 



(7) 



where {4>i a } is a local orbital basis for the system, a subset 
L of which represents localized electrons in the system and 
Y^'ia jp indicates summation with the self term (i,a) = (j, (3) 
omitted. |^o) is the uncorrelated Hartree-like wavefunction 
corresponding to \^g)- The z-factors are renormalization 
weights for the localized part of the one-particle density ma- 
trix 



V n ia (1 n ia) 



(8) 



where p" r , = Pi.rPi,r' I (r' \c-ia | T) | 2 is the probability for a 
transition between two atomic configurations I\ and that 
results in the increase of the occupation of the single particle 
state a at site i by one. The summation is over all configu- 
rations on site ills 18]. For non-localized orbitals Zi a = 1. 
In our present notation, a includes both the orbital and spin 
indices and 



= * 



(9) 



The Zi a are therefore functions of the orbital occupation 
{rii a } and the probabilities pf, r ,. The p^ r , can be expressed 
in terms of the probabilities pi (T) for a local configuration 
r ifl6i The set of local orbitals in L and the set of local 
configurations {I\} with non-zero probabilities are specified 
for the system. 

Under the Gutzwiller approximation, the electron density is 
defined as 



P(r) = X) z ^ z 3^*ia ( r ) <t>iP M (*o 

ia,j/3 
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(10) 



We can also define the localized electron density pi (r) by a 
similar expression, except that the summation is restricted to 
a and j3 in i. 

We will choose E xc [p] to be of the form: 

E xc [p] = Pi Ur + / ^ ^ - W 



Pi (r)i;(r,r')p;(r')d d r(i d r' (11) 



We require our system to be the same as the regular LDA sys- 
tem in the limit when there are no localized electrons. This 
can be achieved if we choose e xc (p) to be the same as in 
LDA. In the limit when all electrons are localized our system 
becomes a multiband Hubbard Hamiltonian. Ur is a sum of 
Slater integrals representing the Coulomb repulsion between 
localized orbitals on the same site in the configuration T. 

The variational degrees of freedom in our system are 
{pi (F)} and l^o}- Since |\I>o) can be expressed as a simple 
product of one particle wavefunctions {V'nk}, it follows that 



*0> = X^"< k (V'nkl&a) (4>jfi\lpnk) ( 12 ) 



where m, k are the usual band indices and /„ t is 1 for oc- 
cupied states and for empty states. The variational parame- 
ters in our calculations are {p; (T)} and {V'nk} with the con- 
straints that {V'nk} are normalized to 1. 

A set of single electron equations can be derived using 
the variational principle by taking the derivatives of the new 
energy functional Eqf2]with respect to {V'nk} and {pi (T)}, 
keeping in mind that the density p (r), the localized density 



3 



pi (r), the exchange-correlation functional e xc (p)and the pa- 
rameters Zi a are defined above. This set of equations can be 
solved self-consistently to give the band structures and total 
energies of the correlated electron system. 

By taking the derivatives with respect to {VVik} we have 

Heftlpnk = Kk^nk (13) 

with effective Hamiltonian 

H cS = H L G + Y, 2e — — ^P iQ (14) 
ia/3 0riia 

where the first term H l G is the Gutzwiller-renormalized oper- 
ator of H where 

H l =f + V ion + V H + fee - Pi (V 1 H + e xc ) Pi (15) 

is the effective mean-field potential with the localized- 
localized electron-interaction contributions subtracted out. 
Vh and Vjj are the mean field Coulomb potential (Hartree 
potential) due to the total and localized charge respectively. 



/'x 



d(p-pi)e xc (p) 
dp 



(16) 



while Pi is the projection operator on the localized subspace 
L, i.e. Pi<f>i a = (f>i a for a 6 L (localized electron orbitals) 
and Pi<pi a = otherwise. 

The second term in Eq[14|adds back the localized-localized 
electron contribution to the effective potential (subtracted 
from the first term) according to the Gutzwiller approxima- 
tion. In the second term, Pi a is the projection operator on <\>i a 
and 



1 /k 



p ia m 



4>nl 



(17) 

The derivatives over the local configuration probabilities 
{jU (T)} yield 



= U r + 2 ^ e i 



d In z ia 



(18) 



This is a set of self-consistency criteria to be satisfied 

byfe (r)}. 

The set of equations in (13) and (18) can be solved iter- 
atively to obtain a self-consistent solution for {pi (T)} and 
{V'nk} and the total energy of the system evaluated according 
to equations (1) to (12). 

The specific derivation of our density functional approach 
is based on the Gutzwiller method. It is, however, possi- 
ble to formally generalize the approach and set up a density 
functional for an arbitrary solvable interacting many-electron 
reference system, independent of the specific aspects of the 
Gutzwiller formalism. To put it in a more general context, we 
discuss these aspects of our approach. A crucial ingredient 



of the density functional formulation by Kohn and Sham is to 
choose the kinetic energy functional T [p] to be the kinetic en- 
ergy of a system of independent electrons in a potential V s that 
yields the ground state density p(r). Our approach differs 
from this key starting point by Kohn and Sham by defining 
T [p] as the kinetic energy of a system of interacting elec- 
trons. We then assume that one can again formulate an effec- 
tive many body problem with Hamiltonian 



H S =T + V S + U S 



(19) 



that yields the ground state density p (r) through an appro- 
priate choice of the single particle potential V s . U s still con- 
tains explicit interactions among the electrons. We then use 
the wave function |^ s ) of this correlated reference system to 
define the kinetic energy functional 



T[p] 



T 



(20) 



This necessitates a new, modified functional for E xc [p] which 
is no longer the only term where correlation effects enter. A 
self-consistent set of density functional equations emerges for 
any choice of |\& s ) that allows for an evaluation of the func- 
tional derivative ST [p] /6p(r). The Gutzwiller wave func- 
tion l^s) = I^g) discussed above or other Jastrow type wave 
functions are examples. The close formal connection to a non- 
interacting electron system (see Eq0]i makes the analysis of 
the kinetic energy feasible and leads to the mapping of the 
many-particle problem onto a set of effective single particle 
Schrodinger equations. 

We can make further progress in our analysis of E xc \p] by 
using the coupling constant integration approach of Ref . II 1 911 . 
We first make a specific choice U s = PUP for the interaction 
term, where P projects onto the configuration space of the 
strongly interacting electrons e.g. those among local 3c?, 4/ 
or 5/ electrons. We then introduce the Hamiltonian 



H x = T + U s + A U - U s 



V(X) 



(21) 



with varying coupling constant A, where the potential V (A) 
equals the nuclear potential Vo for A = 1 and is assumed to 
yield a density p(r) independent of A for A < 1. With the 
help of the Hellmann-Feynman theorem we obtain an explicit 
expression for the exchange correlation functional: 



E xc [p] = (* s 



* s ^+ / d 3 r {pe xc (r) - Pl e XCj i (r)) 
p l {v)v{r,r')p l {Y r )d\d i r' 1 (22) 



that is independent of the specifics of the above Gutzwiller 
approach. The exchange correlation potential e xc j (p) of the 
localized orbitals is found to be 



e xc ,l (r) = \J d3r 'Pi ( r « ( r . O iSl (r, O - 1) . 



(23) 
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where 9 «(,r<) . J ^^^^ is 
the two particle correlation function of localized states with 
b aa i (r, r) = ?/v (r') (r). Comparing this result with 



and 



EqUH we identify £ ir p*(r)[/ r = 
find that we made the approximate choice e XCj ; (p) ~ 
£a;c(p)- Then, the second term in Eq|22] becomes simply 
/ d 3 r (p — pi) e xc (p). This simple choice reproduces cor- 
rectly the limit of the ordinary density functional theory of 
Kohn and Sham in the case without localized orbitals. A sim- 
ilar form was also shown to be very successful in describing 
the non-linear exchange-correlation interactions between core 
and valence charge densities ll2(il . 

An important physical consequence of our new approach is 
its ability to combine two, seemingly distinct, mechanisms for 
screening the Coulomb interaction between electrons. Usu- 
ally, screening is understood as a response of the particle 
density in the vicinity of a charged object and many aspects 
of it are appropriately incorporated in the usual Kohn-Sham 
density functional formalism. On the other hand, in case of 
strong electron-electron correlations, the Coulomb interaction 
can also be screened via a reorganization of the many body 
state, for example in the form of strong band renormalizations, 
amounting to a drastic change in the kinetic energy of the elec- 
trons. These effects are at best poorly described in the usual 
density functional formalism. In our approach, through the 
self-consistent solution for {ipnu} and {pi (T)}, the system 
can respond to the addition of extra terms in the Hamiltonian 
(both in the external potential as well as in Up)- The approach 
has the advantage of combining both aspects of screening in 
a self consistent way without resorting to model parameters 
or model Hamiltonians. This opens a new perspective for the 
first principles description of strong electron correlations in 
complex materials. 

In summary, we have developed a new density functional 
theory incorporating strongly correlated electronic effects into 
the kinetic energy functional within the Gutzwiller approx- 
imation. We show that a set of single particle equations 
can be obtained from functional derivatives of the energy 
with respect to the orbitals included in the non-interacting 
one-particle density matrix with renormalized kinetic and po- 
tential operators. This set of equations can be solved self- 
consistently in a way similar to regular LDA calculations. 
In our scheme, the variational parameters are the single par- 
ticle electron wavefunctions for the localized and delocal- 
ized electrons and the occupancy of the various localized 
configurations at each atom in the unit cell. Like LDA, 
and unlike the many other correlated electron calculations on 
the market, the formulation is from first principles with all 
Coulomb and exchange-correlation interactions determined 



self-consistently. There are no adjustable parameters. We be- 
lieve developments along this approach will be fruitful in ex- 
tending the successful applications of density functional cal- 
culations to new systems with important electron correlations. 

We are grateful to V. Antropov, B. N. Harmon, W. Weber 
and J. Buenemann for useful discussions. Ames Laboratory 
is operated for the U.S. Department of Energy by Iowa State 
University under Contract No. DE-AC02-07CH1 1358. This 
work was supported by the Director for Energy Research, Of- 
fice of Basic Energy Sciences. 



[1] P. Hohenberg and W. Kohn, Phys. Rev. 136 B864 (1964). 
[2] W. Kohn and L. J. Sham, Phys. Rev. 140, Al 133 (1965). 
[3] Quantum Theory of Real Materials, Eds. J. R. Chelikowsky and 

S. G. Louie, (Kluwer Press, 1996). 
[4] VASP (Vienna ab-initio simulation program), G. Kresse and J. 

Hafner, Phys. Rev. B 47, 558 (1993); ibid. 49, 14 251 (1994); 

G. Kresse and J. Furthmuller, Comput. Mat. Sci. 6, 15 (1996); 

G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11 169 (1996). 
[5] V. I. Anisimov, and O. Gunnarsson, Phys. Rev. B 43, 7570 

(1991);V. I. Anisimov, J. Zaanen, and O. K. Anderson, Phys. 

Rev. B 44, 943 (1991); V. I. Anisimov, I. V. Solovyev, and M. 

A. Korotin, Phys. Rev. B 48, 16929 (1993). 
[6] G. Kotliar and D. Vollhardt, Physics Today, March 2004, p. 53 

(2004). 

[7] G. Kotliar, S. Y. Savrasov, K. Haule, V.S. Oudovenko, O. Par- 
collet, and C. A. Marianetti, Rev. Mod. Phys. 78, 865-951 
(2006). 

[8] K. McMahan, Phys. Rev. B 72, 115125 (2005). 
[9] S. Y. Savrasov, G. Kotliar, and E. Abrahams, Nature 410, 793 
(2001). 

[10] L. Petit, A. Svane, Z. Szotek, and W. M. Temmerman, Science 
301, 498 (2003); Eur. Phys. J. B 25, 139-146 (2002). 

[11] N. M. Kudin, G. E. Scuseria, and R. L. Martin, Phys. Rev. Lett 
89, 266402 (2003). 

[12] I. D. Prodan, G. E. Scuseria, and R. L. Martin, Chem. Phys. 3, 
225 (2005). 

[13] M. Gutzwiller, Phys. Rev. Lett. 10, 159 (1963) ; Phys. Rev. 137, 

A1726 (1965). 
[14] D. Vollhardt, Rev. Mod. Phys. 56, 99 (1984). 
[15] G. Kotliar and A. E. Ruckenstein, Phys. Rev. Lett. 57, 1362 

(1986). 

[16] V. Dorin and P. Schlottmann, Phys. Rev. B 47, 5095 (1993). 
[17] J. P. Julien and J. Bouchet, Physica B - Conden. Matter, 359, 
783 (2005). 

[18] J. Biinemann, F. Gebhard, and W. Weber, Foundations of 
Physics, 30, 201 1 (2000); J. Biinemann, W. Weber, and F. Geb- 
hard, Phy Rev. B 57, 6896 (1998). 

[19] J. Harris, and R. Jones, J. Phys. F 4, 1170 (1974). 

[20] S. G. Louie, S. Froyen and M. L. Cohen, Phys. Rev. B 26, 1738 
(1982). 



